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Abstract. Through a linear stability analysis, we investigate the effectiveness of 
a noninvasive feedback control scheme aimed at stabilizing traveling wave solutions 
j^^iKx+iujt q£ ^]-^g one-dimensional complex Ginzburg Landau equation (CGLE) in the 
Benjamin-Feir unstable regime. The feedback control is a generalization of the time- 
delay method of Pyragas J], which was proposed by Lu, Yu and Harrison |2j in the 
setting of nonlinear optics. It involves both spatial shifts, by the wavelength of the 
targeted traveling wave, and a time delay that coincides with the temporal period 
of the traveling wave. We derive a single necessary and sufficient stability criterion 
which determines whether a traveling wave is stable to all perturbation wavenumbers. 
This criterion has the benefit that it determines an optimal value for the time-delay 
feedback parameter. For various coefficients in the CGLE we use this algebraic stability 
criterion to numerically determine stable regions in the {K, p)-parameter plane, where 
p is the feedback parameter associated with the spatial translation. We find that 
the combination of the two feedbacks greatly enlarges the parameter regime where 
stabilization is possible, and that the stability regions take the form of stability tongues 
in the {K, p)-plane. We discuss possible resonance mechanisms that could account for 
the spacing with K of the stability tongues. 
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1. Introduction 



A current challenge to our understanding of pattern-forming systems lies in our ability 
to control the spatio-temporal chaos that many of these systems naturally exhibit. The 
mathematical existence of a plethora of simple spatial or spatio-temporal patterns for 
many nonequilibrium, spatially extended systems is well-established on the basis of 
equivariant bifurcation theory 3 . However, these simple patterns often prove to be 
unstable in a given system, which evolves instead to a state of spatio-temporal chaos. 
Perhaps the simplest and best studied manifestation of spatio-temporal chaos is that 
associated with the one-dimensional complex Ginzburg Landau equation (CGLE) |3], a 
universal amplitude equation that describes spatially-extended systems in the vicinity of 
a Hopf bifurcation. In the so-called Benjamin-Feir unstable regime, the simple solutions 
{e.g. traveling plane waves and spatially-homogeneous oscillations) are all unstable 
to long-wave perturbations. The focus of this paper is a linear stability analysis of 
Benjamin-Feir unstable traveling wave solutions of the CGLE in the presence of a 
feedback control scheme that was originally proposed by Lu, Yu and Harrison j2] in 
the setting of nonlinear optics. The feedback approach is noninvasive, meaning that the 
feedback signal decays to zero once the targeted traveling wave state of the system is 
realized. 

Feedback control methods aimed at stabilizing the unstable periodic orbits 
associated with low- dimensional chaotic attractors have been extensively investigated 
for more than ten years now, and have proved especially effective for nonlinear 
optical systems. The simple approach taken by Ott, Grebogi and Yorke [S] relies on 
applying small perturbations to a system parameter that help maintain the system 
in a neighborhood of the desired unstable periodic orbit. However, this approach, 
which requires an active monitoring of the state of the system so that the feedback 
is appropriately adjusted, can prove impractical in systems that evolve too rapidly. 
Autoadjusting feedback control methods, in which the feedback is based upon current 
and past states of a system, have proven useful for rapidly evolving systems because they 
adjust automatically to rapid changes of the system and require no active monitoring. 
One autoadjusting feedback technique that has attracted considerable attention was 
introduced by Pyragas p. In this approach the feedback is proportional to the 
difference between the current and past states of a system, i.e. the feedback is 
F = ■y{x(t) — x(t — At)), where At is the period of the targeted unstable periodic orbit. 
The method possesses a couple of properties which make it attractive experimentally. 
First, if a state with the desired periodicity is stabilized the feedback term vanishes and 
control is achieved in a reasonably noninvasive manner. Second, this type of feedback 
control may be straightforward to implement in the laboratory when feedback loops are 
practical. The method has been implemented successfully in a variety of experimental 
systems including electronic EI , laser |HI , plasma CHI HI] , and chemical systems 
[T2| ITSt IT^ . A number of modifications of the method of Pyragas have also been 
investigated. For instance, Labate et al. [T3] added a filter to their time delayed feedback 



scheme for a CO2 laser and successfully stabilized periodic behavior; the fi Iter rejected 
one of the characteristic frequencies associated with a quasiperiodic route to chaos in this 
system. In yet another direction Socolar et al. [16^ proposed a method of "extended time 
delay autosynchronization" which incorporates information about the state of the system 
at many earlier times t — nAt (prescribed by positive integers n). This extended time 
delay feedback can successfully stabilize traveling wave solutions of the one- dimensional 
CGLE [T7j, although it fails in two- dimensions [THj . 

For spatially extended pattern forming systems, modifications of the time-delay 
autosynchronization scheme of Pyragas have been proposed which take into account 
not only the temporal periodicity, but also the spatial periodicity of the targeted 
pattern. Perhaps the simplest such scheme is the one proposed in |2], which is the 
one we investigate here. Specifically, Lu, Yu and Harrison used numerical simulations 
of the two-dimensional Maxwell-Bloch equations describing a three level laser system to 
demonstrate that spatio-temporal chaos ("optical turbulence") could be eliminated by 
applying a linear combination of time-delay feedback and an analogous feedback term 
involving a spatial translation, i.e., F = p(i?(x+xo) — -E(x)), where xq is the translation 
vector associated with the feedback. Other schemes for controlling spatio-temporal 
patterns have incorporated spatial filters in the time-delayed feedback. For instance, 
Bleich et al. [19] showed, through linear stability calculations and numerical simulations, 
that traveling wave solutions of a model of a spatially extended semiconductor laser 
could be stabilized using a combination of a Fourier filter and an extended time-delay 
scheme involving multiple time delays. Baba et al. [201 considered problems for which 
a Fourier filter is not appropriate; they incorporate a spatio-temporal filter in their 
time-delay feedback that is derived from a linear stability analysis of the targeted 
unstable periodic orbit in the uncontrolled problem. Other successful approaches to 
feedback control of spatial patterns in nonlinear optics have relied on the Fourier filters 
alone [211 1221 I2ni 1211 I2n|- Numerous investigations of oscillatory chemical patterns 
have incorporated a global time-delay feedback control [211 123 I2H1 1211 • In this case 
the magnitude of the global feedback depends on the spatial average of a chemical 
concentration at an earlier time. Global delayed feedback has been shown to suppress 
spatio-temporal chaos associated with the CGLE jHOl IH!] • 

The goal of our paper is to provide a detailed analysis of the spatial and temporal 
feedback control scheme proposed by Lu, Yu, and Harrison |2I. In particular our goal 
is to derive conditions under which this method can stabilize traveling waves in the 
Benjamin- Feir unstable regime of the CGLE. The analysis is complicated because of 
the time-delay, and because we must determine stability with respect to perturbations 
at all spatial wavenumbers. On the other hand, it is simplified because the exact 
traveling wave solutions of the CGLE (both with and without feedback) take the 
particularly simple plane wave form /^e*^^'"*"*"* (where R and ut are simple functions 
of the wavenumber K). Moreover, as is the case for the CGLE without feedback, the 
plane wave solution enters the linear stability analysis only parametrically through the 
wavenumber K. Specifically, the linear stability analysis reduces to an investigation 



of a pair of complex linear delay-differential equations whose coefficients depend on 
the wavenumber K of the targeted plane wave solution and the wavenumber Q of 
the perturbation but not explicitly on the spatial variable x (nor on time). Our 
main linear stability result is an algebraic condition that determines whether or not 
there exists a time-delay feedback term that will stabilize the traveling wave (for a 
given spatially-shifted feedback term). Interestingly, our analysis identifies an optimal 
time delay feedback to use in determining whether the spatial and temporal feedback 
control method will stabilize the traveling wave. Our algebraic condition for stability, 
which must be checked numerically, reveals that the combination of temporal-delays 
and spatial-shifts in the feedback greatly enlarges the regions where stabilization is 
possible. In some instances we are able to understand the failure of the feedback 
control scheme in terms of simple resonant conditions between the wavenumber K of the 
underlying traveling wave pattern and the wavenumber Q of the unstable perturbation, 
or between the frequency uj of the traveling wave and the frequency associated with 
the perturbation. These resonances can be tuned to some extent using the spatially- 
shifted feedback term and thus are not fundamental obstacles of the type that can be 
present when time-delay feedback alone is used IBHl 01] • (Another approach to 
eliminating fundamental limitations of time-delay feedback in low- dimensional chaotic 
systems has been proposed and investigated by Pyragas |SS]-) We expect our paper 
to contribute to the body of literature in which the success of time-delay feedback 
methods is investigated through a detailed linear stability analysis of the associated 
delay-differential equations EH EHl 1321 IHO] • 

Our paper is organized as follows. In section |21 we describe the feedback to be 
investigated, and review the basic properties of traveling wave solutions of the CGLE 
without feedback. This section also reviews the stability analysis of a related problem, in 
which a time delay feedback term was applied to a system near a Hopf bifurcation |4Uj . 
Section IHl sets up our linear stability analysis for the traveling wave solutions of the 
CGLE in the presence of feedback. We then prove our main stability result which 
allows us to determine the conditions under which the feedback parameters can be 
chosen so that Benjamin-Feir unstable traveling waves are stabilized. In section |31 we 
apply our stability criterion to numerically determine the parameter regions in which 
feedback control is possible, and discuss possible cause for failure. Finally, in sectional 
we summarize our results and suggest possible extensions of our analysis. 

2. Background and Problem Formulation 

The 1-D complex Ginzburg Landau equation, an amplitude equation describing a 
spatially extended system near the onset of a Hopf bifurcation, is given by 




(1) 



where A{x, t) is a complex amplitude, and hi and 63 are real parameters. We consider 
traveling wave solutions of (Q) of the form 

Atw = i?e*^"+'"*, (2) 

where the amplitude R and frequency u are determined by the wavenumber K through 
the relations 

- ^i^. (3) 

00 = R^- biK^. (4) 

We focus on the case of a supercritical Hopf bifurcation so that 63 > 0; thus li^l < 1 
follows from Q. 

In the Benjamin-Feir unstable regime {bi > 63 > 0), all plane wave solutions of the 
form (j2I) are unstable. In the Benjamin-Feir stable regime (61 < 63, 63 > 0), traveling 
waves with wavenumbers K satisfying < K"^ < K^^^ = (63 — bi)/{3bs — 61 + 2/63) 
are linearly stable to longwave perturbations 021 HHj- Figure Q summarizes the 
behavior of typical solutions of in the {bi, 63)-parameter plane. The Benjamin-Feir 
unstable regime is divided into three different regions, an amplitude turbulent regime, 
a phase turbulent regime, and a bichaotic regime. The amplitude turbulent regime 
is characterized by large fluctuations in amplitude and phase defects. In the phase 
turbulent regime amplitude fluctuations are much smaller and no defects are observed. 
In the bichaotic regime either amplitude or phase turbulence occurs depending upon 
initial conditions |3]. 

We consider how feedback of the type proposed in j2] affects the hnear stability 
of the plane wave solutions of the CGLE in the Benjamin-Feir unstable regime. In 
particular, we consider 

dtA = A + {l + tbi)dlA - (63 - i)WA + F (5) 

where the feedback term, F, is a linear combination of spatially translated and time- 
delay feedback terms. Specifically, 

F = p[A{x + Ax, t) - A{x, t)] + 7[A(x, t) - A{x, t - At)], (6) 

where 7 and p are real parameters and At > 0. Note that if a solution A{x,t) is 
spatially periodic with period L = Ax/n, n & Z and temporally periodic with period 
T = At/m, m & Z , the feedback term vanishes. Thus such a space and time periodic 
state is also a solution of the CGLE without feedback. We will restrict our linear 
stability analysis to traveling plane wave solutions Q which have the same spatial and 
temporal scales as the feedback, i. e. we assume 

Ax = 2n/K , At = 2n/\u;\ . (7) 

Note that this choice of spatial feedback scale assumes that the targeted traveling wave 
state is not too close to bandcenter {K = 0). The temporal feedback assumes a; 7^ 0, or 
equivalently 7^ 1/(1 + 6163). 
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Figure 1. Phase diagram for tlie 1-D complex Ginzburg Landau equation |^ in 
the (63, 6i)-parameter plane. All plane wave solutions Q are unstable above the 
Benjamin-Feir (BF) line bi =63. This figure is reproduced (with permission of the 
author) from |44|. 




Figure 2. Example of critical curves (solid lines) in the {a, G)~parameter plane of ((SJl, 
given by Q, in the nondegenerate case PAt ^ mr for any integer n. (Here we used 
Ai = 1 and /? — — 47r/5.) The dashed lines are the asymptotes G — zta. The numbers 
in the various regions indicate how many solutions of the characteristic equation @ 
have positive real part. The stable region is to the left of all critical curves [10] . 
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Figure 3. Examples of critical curves associated witli © in the degenerate cases for 
which PAt = nir for some integer n. (a) At — 1, (3 ~ — tt (representative of n odd). 
(b) At — 1, f3 ^ —2tt (representative of n even). 



(a) 



Our linear stability analysis builds on some results of Reddy et al. jlO], who 
considered the effects of a time-delayed linear feedback term added to a Hopf bifurcation 
normal form. Of relevance to our analysis are their results on the linear stability of the 
origin, which is determined by an analysis of 

i(t) = {a + i(3)z{t) - Gz{t - At) , (8) 

where a, (3, G, At are real parameters with At > 0. The characteristic equation 
associated with solutions of the form z = e^*/^* is 

A = (a + i(3)At - GAte-^ . (9) 

There are an infinite number of solutions A to this transcendental equation; if all of 
them satisfy Re(A)<0, then 2; = is stable. If any solution A has positive real part, then 
2; = is unstable. (See, for example, Diekmann et al. [45J or Driver [IH] for background 
on delay equations.) 

A solution A of the characteristic equation is purely imaginary on the "critical 
curves" in the (a, G)-parameter plane. These curves are significant because they 
represent boundaries in the [a, G)-plane across which the number of solutions with 
Re(A) > changes. Figure |21 shows an example of critical curves in the (a, G')-plane, 
associated with a "nondegenerate case" for which f3At 7^ nvr for any n & Z, with 
examples of critical curves for degenerate cases provided in Figure El A parameterization 
of the critical curves is obtained by letting A = ii/, z/ G -R, in (jU)), in which case we find 

(z/ - (3 At) cos(z/) 



a 



At sin(u) 



G^ 'r^? - (10) 

Atsm(i^) ^ ^ 

Here each critical curve is associated with an interval u G (mvr, (m+l)7r) for some integer 
m. In the nondegenerate situation, Reddy et al. |4Q|, extending a result in show 



that the stabihty boundary is determined by the left most critical curve in the {a, G)- 
plane, as indicated in Figure |21 This stability boundary corresponds to the critical 
curve for the i/-interval containing /3At; it passes through the origin and approaches 
the asymptotes G = ±a as a — — oo. In the degenerate cases (Figure 01), defined by 
(3 At = riTT for some integer n, there is an additional critical curve coincident with one of 
the asymptotes; specifically if n is even (odd) then on the critical curve a = G {a = —G) 
the characteristic equation Q is solved by A = zz/ = il3At{= inir). Finally, we show 
that the following properties hold for the critical curves in the nondegenerate case. 

Lemma 1. Consider G < and /3At ^ nn for any n E Z in (0)- (1) The stability 
boundary in the {a, G) -parameter plane is the greatest distance from the asymptote 
a = G for G = (2) All other critical curves (to the right of the stability boundary 

and for G < 0) lie in the region G < 

Proof. (1) At a point where the distance between the stability boundary and the 
asymptote a = G is maximized 

da dG _ d'^G d'^a ^ ^ 

dv du ' dv"^ du'^ ' 
where cx{v) and G{v) are the parametric equations for the critical curves ()10|) . It follows 
by direct calculation that these conditions are met when = z/*, where v* is defined 
imphcitly by the equation 

sin(z/*) = 13 At - V* . (12) 

From (jmi), we see that this implies G = (It follows from property (2) of this 

lemma that v* defines a point on the stability boundary rather than another critical 
curve.) (2) The remaining critical curves in the region G < reach their maximum in 
the (a, (j)-plane at points for which ^ = (while ^ 7^ 0), where 

dG sin(z/) — (z/ — (3 At) cos(z/) , , 

dp ~ Atsin^{v) 

Substituting 

sin(z/) = (z/ - [3 At) cos(z/) (14) 
into (Uni) yields 

Gmax = -T- -J-T (15) 

At cos(z/) 

for the maximum value of G, where z/ satisfies ()14j). For the critical curves confined to 

G < 0, -1 < COs(z/) < in dH)) so Gmax < □ 



3. Linear Stability Analysis 



3.1. Preliminaries 



In order to determine whether the feedback (0171) can stabilize a traveling wave 
solution (j2I) of (jni), we analyze the effects of small amplitude perturbations on the 
solution by letting 

A = i?e*^^+*^*(l + a+(t)e*Q^' + a.{t)e-'^'^). (16) 

Here a_|_ and a_ are the amplitudes of small perturbations with wavenumbers K ± Q. 
(We may, without loss of generality, assume that Q > 0.) Substituting (fTBj) into (0) and 
linearizing, yields the following system of delay equations: 





(17) 



where 



J 



7-Q(2ir + Q)(l+26i) 



-{h-i)B' 



+ Q{2K ~ Q){1 - th) 
-{b3 + i)R^ + p[e'^^'' -1] 

Diagonalizing yields two decoupled linear delay equations. 



d . , 



jkCkit) ~ lCk{t - At) , A; = 1,2, (19) 

where ji and j2 are the complex eigenvalues of J defined such that -Re(ji) > Re{j2). 

Note that the temporal and spatial feedback terms, proportional to 7 and p, simply 
shift the eigenvalues of J by the diagonal factor 7 + p[e*'^'^^ — 1] . In particular, each jk is 
of the form = 'J+jk, where jk is independent of 7 and depends on 61, 63, p, K, and Q. 
In the following we denote the real and imaginary parts of jk by jkr and jki, respectively. 
Since jir plays a central role in our analysis, we simplify notation by replacing it with 
/ and suppressing its dependence on the system and solution parameters bi,b3,p,K, 
viewing it as a function of the perturbation parameter Q. Specifically, we let 

f{Q) = jZiQ; 61, 63, p, K) = j,,{Q- 61, 63, p, ir, 7) - 7 • (20) 

We note that f{Q) and jii{Q) are continuous functions of Q, and that f{Q) satisfies 

/(O) = 0, 

hm f{Q) ^ - 00 . (21) 

Q—*oo 

Hence f{Q) has an absolute maximum fmax > for some Q > 0. The significance 
of f{Q) is that it represents the growth rate of perturbations of wavenumber Q in the 
absence of the time-delay feedback {i.e. for 7 = 0). 

The characteristic equations associated with solutions Ck{t) = e'^'=*/^* of fT^ are 



Afc = At{jk 



76 



A; = 1,2 



(22) 



We obtain the critical curves associated with our hnear stabihty problem (|19|) by 
substituting = iuf^ into ((221) • Then, for instance, we find the following parametric 
representation of critical curves in the (ji^, 7)-planes: 

. _ (z/i -juAt) cos(z/i) 
At sm(z/ij 

_ (^^IH^. (24) 

Similar equations apply for the critical curves in the (j2r; 7)^parameter plane. Using 
jlr = 7 + f{Q) (which, recall, defines f{Q)) we determine that on the critical curves in 
the (jlr, 7)-plane Q and vi are related through 

f{Q)M = w{ur,Q), (25) 

where 

I _ (//i-Ji.(g)At)(cos(z/i)-l) 

w[vi,Q) = — , (26) 

sm(z/i) 

which is a result we use later in our analysis. 

The parametric equations for the critical curves ()23H24j) in the (jir? 7)^plane are 
similar to those obtained by Reddy et al. jlO], given by (fTn|l . However since both jir 
and iu depend on the perturbation wavenumber Q the results in jlHI, summarized by 
Figures El and El do not apply directly. Nonetheless some useful results concerning 
the stability boundary do carry over. For instance, it follows from ()23|) and (j^^ that 
Jlr = 7cos(z/) so all critical curves must lie in the region of the (jir,7)-plane where 
|jir| < |7|- In other words, if we divide the (ji,., 7)-plane into four quadrants separated 
by 7 = ±jir the critical curves lie inside (or possibly on the boundaries) of the upper 
and lower quadrants (as is also the case in Figures El and E]) . The stability assignments 
are necessarily the same at all points inside the left and right quadrants of the parameter 
plane. In particular, the left quadrant is contained in the stable region and the right 
quadrant is contained in the unstable region for all ju. Because the stability boundary 
is contained in the upper and lower quadrants both stable and unstable regions are 
located there, as signified by the question marks in Figure IU For fixed values of 6i, 63, 
and p, the quadrants of the (jir,7)-plane may be mapped to the ((5,7)-plane. In 
particular, the asymptotes 7 = ±jir- are mapped to the curves defined by f{Q) = and 
7 = —f{Q)/2; an example is shown in Figure 01 

We note that due to the translation symmetry of (0), which acts nontrivially on the 
traveling wave solution there is always a neutral mode associated with perturbations 
at Q = 0. (It is this symmetry that forces /(O) = 0.) Thus the line jir = 7 is a 
(degenerate) critical curve; as it is a consequence of translation symmetry, it does not 
necessarily herald an instability. 

We now make use of the elementary (in) stability results summarized by Fig. Efa) 
to prove Lemmas EHSl which allow us to restrict any further linear stability analysis to 
the k = 1 delay equation of (fTI?|) for Q > 0, and for feedback parameters 7 < and 




Figure 4. Mapping of tlie guaranteed stable (S) and unstable (U) regions bounded 
by 7 = ±jir from (a) the (ji^, 7)-plane to (b) the (Q, 7)-plane. The parameter values 
are bi = 5, 63 = 1.2, K = 0.18, and p 0.1. 

p > 0. In particular, we find that the restriction 7 < 0, p > is sufficient to ensure 
that the feedback does not further destabilize the traveling waves in the Benjamin- Feir 
unstable regime, i.e. the feedback does not introduce new instabilities in this case. 

Lemma 2. (a) If f{Q) < for all Q > 0, where f{Q) is defined by the relation / f^) . 
then the traveling wave is stabilized for any 7 < 0. In particular, the wave is stabilized 
using spatial feedback alone in ^ (i.e. for 7 = Oj, and there are no instabilities 
associated with (5 = when 7 < 0. (b) If f{Q) > for any Q then the traveling 
wave cannot be stabilized using 7 > 0. 

Proof, (a) If f{Q) < 0, then j2r < jir < 7 for all Q > since jir = 7 + f{Q)- Hence, for 
7 < 0, the system is in the guaranteed stable regime of Fig. I^a) for all perturbations 
with Q > 0. It remains to eliminate the possibility of instabilities at Q = 0, where 
J2r(0) < jir(O) = 7. For this we consider the solutions = fik + i^^k of the characteristic 
equation (j^ . and determine that the growth rates p^ satisfy 

Pfc = At(£r + 7(1 - e-^'^- cos(i/fc))) . (27) 

There are no p^ > solutions of (P7j) when 7 < and jkr < 0. This is proved by 
assuming there is a solution with p^ > 0, which would imply that the left- and right- 
hand-sides of ipTj) do not have the same sign, thereby leading to a contradiction. In 
fact, the only solution with p^ = is the one forced by translation symmetry. Hence 
there can be no linear instabilities associated with Q = when 7 < 0. 

(b) If 7 > and f{Q) > for some Q, then the system is in the guaranteed 
unstable regime of Fig. EJ defined by 7 > —^f{Q), f{Q) > 0, for perturbations with 
that wavenumber. If there is no temporal feedback (7 = 0) and f{Q) > for some Q, 
then instability for these Q values follows directly from the ci{t) equation of (fT^ . which 
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^ = {f{Q)+^Jl^)cl for 7 = 0. (28) 

□ 

Lemma 121 motivates our focus on 7 < in the remainder of the paper: if the time- 
delay feedback is to be effective in stabihzing plane waves, then the value of 7 used 
in (jH)) should be negative. We similarly focus on p > 0. This is because the spatial 
feedback has the effect of shifting each jkr in (j27jl by p[cos{QAx) — 1]. Thus there is 
no effect of the spatial feedback on perturbations with wavenumber Q = 2n7v/Ax for 
some n & Z. And, for all other wavenumbers, there is a shift towards the stable region 
of Figure Efa) if p > and towards the unstable region if p < 0. In particular, in 
the absence of time-delay feedback, if p < the spatial feedback exacerbates longwave 
instabilities in the Benjamin-Feir unstable regime. In the remainder of the paper we 
consider only p > 0. In this explained by the following lemma, we may ignore 

the C2-equation of ()19|) in the remainder of our linear stability analysis. 

Lemma 3. // 7 < and p > in / lj.9)) then all solutions of the C2-equation of 
decay exponentially. 

Proof. The real parts of the eigenvalues jk of the matrix J ()18|1 . can be written in the 
form 

ji^ = a + 7 + \p\ , jar = a + 7 - IpI (29) 
since they are solutions of a quadratic. From p8|) it follows that 

« = -(Q' + &3i?' + p[l-cos(gAx)]), (30) 
so a < for p > 0. Thus 

J2r-l = a-\p\<0 (31) 

for all values of Q and p > 0. From Fig. |H^a) (with j^r replacing ji,.) we see that the 
parameters are in the guaranteed stable region when 7 < 0. Finally, if 7 = 0, then the 
C2-equation reduces to 

= j2C2 for 7 = 0, (32) 

dt 

so all solutions decay for 7 = since j2r = a — \p\ < 0. □ 



3.2. Main Stability Results 



Lemma|2lstates that if /(Q), defined by pO|l . is negative for all Q > 0, then the traveling 
wave may be stabilized using spatial feedback alone (7 = 0). This section addresses the 
more difficult case where the chosen spatial feedback fails to stabilize the traveling wave, 
i.e., we consider the situation where f{Q) > for some range of Q values and investigate 
whether the traveling waves may then be stabilized by the addition of an appropriate 
time-delay feedback with 7 < 0. (From Lemma El we know that the wave cannot be 
stabilized with 7 > 0.) We focus on the case where the spatial feedback parameter p 
is fixed at some non-negative value so that Lemma El applies. We show that it is then 
sufficient to consider just 7 = — 1/At in order to determine whether or not the traveling 
waves can be stabilized by including the time-delay feedback in 

Our analysis relies on the observation that not all regions of the {jir, 7)-parameter 
plane are accessible for fixed 61, 63, K, p and Q > 0. In particular, since jir = 7 + f{Q) 
and f{Q) has an absolute maximum (denoted by fmax), parameter values to the right 
of the line 

7 = Jir - fmax ("^2) 

are never achieved for any value of Q (see Figure E}- We refer to the boundary (jHSl) 
between accessible and inaccessible regions of the (ji^, 7)^parameter plane as the 
"existence line", meaning that there exist values of Q that place ji^ and 7 on or 
to the left of this line, but never to the right of it. Our stability analysis focuses 
on determining whether there is a (negative) 7-value such that the traveling wave is 
stable for all values of jir between —00 and the existence line, i.e. for all values of Q. 
Fig. El presents, schematically, projections onto the (ji^, 7)-plane of the points that lie 
on the stability boundary in the three-dimensional (jir(Q), 7, jii(Q)) parameter space. 
Figure ini^a) typifies the case where it is possible to stabilize the traveling wave over the 
range of 7-values for which the stability boundary leaves the existence region. On the 
other hand, in Fig. Efb) the stability boundary lies entirely to the left of the existence 
line for all values of 7 so it is not possible to stabilize the traveling wave in this case. 

We now state our main linear stability result; the proof, which relies on Lemma |H1 
below, is deferred to the end of this section. 

Theorem 1. The traveling wave solution ^ of the CGLE ^ with feedback may 
be (linearly) stabilized, for a fixed given value of p > 0, if and only if the following 
equations fail to have a real solution {u, Q), Q > 0: 

u + sm{u) = At juiQ) (34) 

1 - cos(z/) = At f{Q) . (35) 

// a solution fails to exist then the traveling wave can be stabilized using 7 = —1/At. 

Remark 1: Note that any solution of (jH^j) that exists must lie in a closed interval 
of Q values defined by the requirement f{Q)At E [0, 2]. 



0.2 0.4 0.6 0.8 1 ■" 

(a) (b) 

Figure 5. (a) Plot of f{Q) for the parameters of FigureQl (b) Corresponding regions 
of the (jlr, 7)-parameter plane (where jir = 7 + f{Q))\ jir values to the right of the 
existence line are inaccessible. 




Figure 6. (a) Schematic of stable region for case that the stability boundary leaves 
the existence region. Range of 7 values for which the traveling wave is stable against 
perturbations of all wavenumber Q is indicated, (b) Schematic of the unstable case 
where the stability boundary lies to the left of the existence line for all 7. 

Remark 2: Equations (j3lH35|l are equivalent to the equations for the critical 
curves (|2;-ill24jl for 7 = -1/At and Ui = u. 

The following lemma addresses special degenerate situations for which stabilization 
of the traveling wave is impossible. 

Lemma 4. Let Q = Q2n > be defined by the condition that jii{Q2n)At = 2?™, n E Z. 
The traveling wave solution (0j of the CGLE (Op with feedback is linearly unstable 
for all values of'y if for any n (a) f{Q2n) > 0, or (b) f{Q2n) = and ^{Q2n) 7^ 0. (c) 
If f{Q2n) = ^{Q2n) = 0, then the traveling wave is either linearly unstable or neutrally 
stable; it is not linearly stable for any value of 

Proof, (a) In the degenerate situation, for the value ofQ = Q2n, where jii{Q2n)^t = 27m 
{n G Z), the hne jir = 7 represents a critical curve for the delay equation (fTTI|l (with 



k = 1). (Here we are viewing ji,. as an independent parameter, with ju held fixed; see 
Fig. Efb) with a = jir and G = 7.) Thus the travehng wave solution is unstable to 
perturbations with wavenumber Q = Q2n if 3ir{Q2n) > 7- Since f{Q) = jir{Q) — 7, this 
inequality is met whenever f{Q2n) > 0. 

(b) If f{Q2n) = 0, then jir{Q2n) = 7 and the solution is neutrally stable to 
perturbations of wavenumber Q = Q2n for 7 < 0. (It is unstable if 7 > since it 
lies in the guaranteed unstable regime of Fig. IH^a)) To determine the stability of the 
traveling wave for 7 < we examine the movement, with Q, of the solutions Ai = fii+iui 
of the characteristic equation ()22j). From the real and imaginary parts of ()22|) we have 

/ii = f{Q)At + 7At - 7Ate~^i cos(i/i), (36) 

z/i = MQ)At + ^Ate~^' sin(i/i), (37) 

which has solution /ii = 0, ui = jii{Q2n)At = 27in for Q = Q2n- To determine the 
movement of /ii with Q near Q2n, we compute ^{Q2n)- We find 

Since j^ki > for 7 < 0, the eigenvalue Ai moves into the unstable region /xi > with 
increasing (decreasing) Q if ^{Q2n) > (< 0). 

(c) If f{Q2n) = ^{Q2n) = 0, then the traveling wave may only be neutrally stable 
for 7 < 0. This latter claim follows since the critical eigenvalue Xi{Q2n) = 2vm does not 
cross the imaginary axis with finite speed as Q is varied near Q2n, i-e. ^{Q2n) = 0. 
This tangency implies that the traveling wave cannot be linearly stable in this case - it 
is either unstable or "at best" neutrally stable. □ 

A consequence of this lemma is that the temporal feedback cannot be used to 
stabilize homogeneous oscillations. (Note that for K = 0, spatial feedback is not 
applicable and only the temporal feedback is relevant.) Specifically, for K = p = 0, 
the matrix J, given by (fTHj) . has eigenvalues that are purely real for Q sufficiently small 
and hence we have jii{Q) = with f{Q) > for sufficiently longwave perturbations 
associated with the Benjamin-Feir unstable regime. It follows immediately from 
Lemma |3] (with n = 0), that the feedback fails to stabilize the uniform oscillatory 
mode associated with K = 0. In this context, we note that Beta et al. jSH] have 
given a possible physical explanation for the failure of time-delayed global feedback to 
stabilize a uniform oscillatory mode in the case of diffusion-induced chemical turbulence. 
Harrington and Socolar ^H] have a related result that shows that time-delay feedback 
alone cannot stabilize traveling wave solutions of the two-dimensional CGLE. 

We may further interpret Lemma 0] in light of analogous results for the spatial 
feedback. Specifically, we recall that the effect of the spatial feedback on the eigenvalues 
of the matrix J{Q), given by p8|l . is to shift them by the diagonal entry p(e*'^'^^' — 1). 
Thus if there is an unstable wavenumber Q = Qn for which Q„Ax = 2mT for some 
integer n, then the spatial feedback can have no stabilizing effect on that perturbation 



wavenumber. In this case, the perturbation wavenumber Q„ is "resonant" with the 
underlying wavenumber K of the travehng wave [Qn = nK), and it does not feel the 
influence of the feedback. Lemma ^ leads to an analogous result when the frequency 
associated with the perturbation at Q2n, namely iu{Q2n)i is in resonance with the 
temporal delay At, i.e. jii{Q2n)^t happens to also be an even multiple of vr. In 
that case, the temporal feedback term also has no effect on the perturbation, i.e. 
it cannot suppress the instability associated with it. Just et al. refer to such 
situations where the method of Pyragas fails as "torsion- free" . These observations 
about special (nongeneric) spatial and temporal resonant cases are summarized by the 
following corollary to Lemma IH 

Corollary 1. //, in the absence of any feedback ('y = p = 0), there is an unstable 
wavenumber Q = Q for which QAx = 2m7r and jii{Q)At = 2mT for some integers m 
and n, then the traveling wave cannot be stabilized using the feedback (0j for any value 
of p or 7. 

Remark: It follows from ((Zj) that the resonance conditions QAx = 2rmT and 
jii{Q)At = 2n7T are equivalent to the conditions Q = mK a.nd jii{Q) = nu, respectively. 

It now remains to consider whether stabilization is possible when jii[Q)At 7^ 27m 
for any n ^ Z over all intervals of Q for which f{Q) > 0. It is consideration of 
this situation that leads to the stabilization criterion of Theorem ^ In particular, we 
aim to determine a necessary and sufficient condition for the stability boundary in the 
(iir, 7)^plane to be parameterized by Q, such that it extends continuously from the 
origin to 7 — * — 00 along the asymptote 7 = jir, as in Fig. IH^a) Hence we focus on the 
mapping ()25fOH|l between Q and Ui along the critical curves. 

We first describe some properties of the function w(z/i, Q), given by (j^ . for fixed 
Q. In particular, we are interested in this function over the interval of z/i-values that 
parameterizes the stability boundary in the (ji^, 7)-plane in the case where jir is allowed 
to vary independently of ju, which is held fixed with juAt 7^ 2n7r in ()23ll24j) . Our first 
lemma treats the case where juAt 7^ nvr for any n ^ Z, while the second lemma 
describes the degenerate case juAt = {2n + l)7r for some n & Z. Each lemma consists 
of two parts: the first part merely defines the interval of z/i-values that parameterizes 
the stability boundary in each (ji^,, 7)-plane {ju fixed) and the second part describes 
w(z/i, Q) on that interval. The first part of each lemma follows directly from the analysis 
of Reddy et al. [40j, described in Section E| 

Lemma 5. Consider •y < and let jir in ^EBl(^ vary independently of ju, which is held 
fixed with juAt 7^ nvr for any n a Z . (1) The stability boundary in the {ju-i'y) -plane 
is associated with an interval 1^ = {vai^h) of vi-values. The endpoints, z/^ and Uf,, are 
ju^t and 27m, where n & Z is chosen so that |jijAt — 27{n\ < tt. (2) For vi G ly the 
function w{vi,Q) (with Q fixed) is positive with w{ui,Q) ^ as z/i approaches z/^ and 
Vh', w{vi,Q) has a unique maximum value Wmax{Q) (it Vi = v* , where 

sin(z/*) = ju{Q)At - u* , WmaxiQ) = 1 - cos(z/*) . (39) 
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Figure 7. Plots of w{i'i,Q) vs. vi, where w is defined by (|26|l . for various (fixed) 
values of Q such that j'l^At g (0, tt]. 

Proof. From Reddy et al. |40j, we have that the stabihty boundary (for 7 < 0) extends 
between the origin of the (ji^, 7)-plane and the asymptote 7 = jir as 7 —00. The 
boundary is associated with an open interval /j^; the endpoints of the interval correspond 
to jijAt (at the origin) and the nearest value of 2?™ (as 7 — > —00). In particular, the 
following two possibilities follow from the parameterization ()23II24|) : (a) there exists an 
n E Z such that juAt — 27m < vr in which case 1^ = (27m, jijAt), or (b) there exists 
an n E Z such that 27m — juAt < 7r in which case 1,^ = (jijAt,27m). In both cases 
it is straightforward to show that w{v,Q), given by (|26|) . is positive on the interval, 
approaching at the endpoints 27m and juAt. 

For fixed Q, we find the maximum of w{vi, Q) on I^, by examining 

dw{ui,Q) ^ {juiQ)At - z/i - sin(z/i)) ^^^^ 
dui (1 + cos(z/i)) 

For Ui E lu, the denominator on the right- hand- side of (pn|) is positive and the numerator 
is a decreasing function of z/i which goes through zero at z/i = z/*, where u* is defined 
implicitly by The solution i/* of ^ = is unique and corresponds to a maximum 

of w on the interval 1^. In particular, it follows from ()40|) that w(z/i, Q) is increasing for 
^1 G (z/q, u*) and decreasing for i^i G (i^*, z/ft). (See FigureElfor sample plots of w{h'i, Q) 
on 4). □ 

Lemma 6. Consider 'y < and let jir in vary independently of ju, which is held 

fixed with juAt = (2n + 1)tt for some n E Z . (1) The stability boundary in the (jir,^)- 
plane consists of the line segment 7 = —jir for 7 G (0, —-^] and the parameterized 
curve \2'J\24^ obtained (equivalently) with either the interval = (2n7r, (2ri + l)7r) or 
/+ = ((2?T,-|- l)7r, 2(?T, + l)7r). (2) For z/i G u^(z/i, Q) is a positive function, with w — * 
as Vi — s> 2m7r (m = n, m = (n + 1) for , , respectively), and w approaches its 
maximum value Wmax = 2 as z/i {2n + l)7r. The function w{vi,Q) on I:j is mapped 



to w(z/i, Q) on /+ under the reflection v 2{2n + l)7r — u. It is an increasing function 
on I~ and hence a decreasing function on . 

Proof. (1) In the degenerate case juAt = {2n + l)7r, the hne 7 = —jir is a critical 
curve which forms the stabihty boundary for 7 > — see Fig. |3fa). This part of the 
stabihty boundary meets another critical curve at jir = —7 = this follows from 
considering ()23II24|) in the limit ui {2n + 1)tt. The latter critical curve extends to 
Jir; 7 —00 along 7 = jir in the limit Ui 2mT,2{n + 1)tt. That the z/i-intervals 
I~ and Jjl map out the same curve in the (jir, 7)-plane follows from the reflection 
symmetry of (j23ll24p under u —>■ 2jijAt — u for juAt = {2n + l)7r. This symmetry is 
also necessarily manifest in w{vi,Q). 

(2) The proof of the other claims about wiyi, Q) on If for juAt = {2n + l)7r are 
straightforward. We note that the monotonicity of w{vi, Q) on is determined by (001) • 
In particular, for jijAt = (2r?, + l)7r, there is a unique zero of 1^ at z/i = (2n + l)7r, which 
corresponds to a maximum of w{vi, Q) on the interval {2mT, 2{n + l)vr). Thus Q) 
is an increasing function on /~, while it is decreasing on the interval Jj". □ 

We may now appreciate the relationship between ui and Q along the stability 
boundary by considering the following procedure: for each Q for which f{Q) > 
calculate jii{Q) and f{Q) and then plot the left and right hand sides of ()25|1 as functions 
of Ui over the appropriate interval I^, to determine which (if any) ui values correspond 
to each value of Q. The left-hand-side of ()25j) yields a horizontal line with positive 
intercept. If jii{Q)At 7^ nn, then by Lemma the right-hand-side of (j23), given 
by (j2ni); is a positive function with a unique maximum (and no other critical points) 
over Iiy] see Fig. [3 for an example. Hence, in this case, there are either 0, 1, or 2 values 
of Ui associated with each Q. In this nondegenerate case we characterize the number of 
intersections of the curves associated with the left- and right-hand-sides of over 
the interval I^ by determining whether the following distance function is positive (no 
intersections), zero (1 intersection) or negative (two intersections): 

D{Q) = f{Q)At-w^^aAQ) , (41) 

where WmaxiQ) is given by (jSHI)- We note that D{Q) is continuous. This follows from 
the continuity of f{Q) and jii{Q)- In particular, since ju{Q) is continuous, the unique 
solution u* of (jH^ is continuous in Q, thereby ensuring the continuity of WmaxiQ)- 

For the degenerate cases prescribed hj Q = Qm, where Qm is defined by 
jii{Qm)At = mn, we define 

wiui = mvr, Qm) = lim w{ui, Qm) = I - (-1)" • (42) 

This corresponds to WmaxiQm), which is 2 for odd m and for even m. We now 
address the case that m is odd, i.e. m = 2n + 1 for some n E Z, in which case 
'w{i'i, Q2n+i) is a positive function on the interval (2n7r, 2{n + l)vr), and it is symmetric 
about ui = {2n + l)7r. Thus, for Q = Q2n+i defined by jii{Q2n+i)At = {2n + l)7r, 
it follows from Lemma El that there are either no solutions or one solution of (1^K|) for 



z/i G (2n7r, (2n + l)7r] or, equivalently, for z/i G [(2n + l)7r, (2ri + 2)n). We recall that 
in this degenerate case there is always an additional point on the stability boundary 
for 7 = — /((52n+i)/2 provided < /((52n+i) At < 2; this point lies at the intersection 
of the lines jir = —7 and jir = 7 + f{Q2n+i) (see part (1) of Lemma IH)). Moreover, 
the condition < /(Q2n+i)At < 2 is equivalent to /(Q2n+i)At < Wmax{Q2n+i) since 
Wmax = 2 in this case. Thus the interpretation of the sign of D{Q), defined by (jUJ, 
is the same in this degenerate CclSG clS it is in the nondegenerate case. Specifically, if 
D{Q2n+i) < then there are two points on the stability boundary associated with this 
value of Q2n+i] if D{Q2n+i) = 0, then there is a single point on the stability boundary 
(given by (jir, 7) = (1/At, — 1/At)); and if D{Q2n+i) > 0, then there are no points on 
the stability boundary associated with this value of Q2n+i- 

The next lemma addresses the continuity of particular critical curves that are 
parameterized by Q. It is specifically concerned with the case that jii{Q)At = (2n + l)7r 
(n G Z) at some point along the curve. At such points the relevant interval I^^ of i>i- 
values may "jump" {e.g. from /~ to Jj" in the case that Q = Q2n+i)- However, this 
does not lead to a corresponding discontinuity in the critical curve parameterized by Q, 
as we now show. 

Lemma 7. Consider a critical curve in the {jir, j) -plane, parameterized by Q, which 
extends from the origin to '~f ^ —00 along the asymptote 7 = jir- The critical curve 
is continuous at any points along it that are parameterized by Q = Q2n+i, where 
jii{Q2n+i)At = {2n + l)7r for some integer n, and f{Q2n+i)At G (0,2). 

Proof. Let Q2n+i be in the set of Q values that parameterizes the critical curve in 
the (jij., 7)-plane, where Q2n+i is defined by the condition jii{Q2n+i)At = {2n + l)7r 
for some integer n. Provided f{Q2n+i)At G (0,2), then there are two points in the 
(jir, 7)-plane, denoted Pi and P2, associated with Q2n+i- (See part (1) of Lemma El) 
Specifically, let Pi be the point that is associated with the unique solution of ()25ll2(')j) on 
/~ = (2n7r, (2n + l)7r) (or, equivalently, on J+ = ((2n + l)7r, (2n + 2)7r), and let P2 be the 
point that lies at the intersection of the lines jir = —7 and jir = /(<52n+i) + 7- Since 
jii{Q)At varies continuously with Q, then for values of Q in a neighborhood of Q2n+i 
we expect the interval of i^i-values associated with the mapping ()25II26|) to be as 
described in LemmaEl For example, if ^{Q2n+i) > 0, then jii{Q)At increases through 
{2n + 1)tt as Q increases through Q2n+i and the relevant interval of i^i-values switches 
from {2mi, jii{Q)At) to {jii{Q)At,2{n + l)7r) as Q increases through Q2n+i- In any 
case, for Q sufficiently close to Q2n+i so that f{Q)At remains smaller than Wmax{Q)^ 
there are two solutions vi of ()25H26p on the appropriate /,^-interval. These values of vi 
determine a pair of points on the critical curve in the (jir, 7)-plane via the parametric 
equations ()23II24|) . We now show that these points converge to Pi and P2 as Q — Q2n+i- 
Specifically, one of the pair of solutions is associated with a i^i-value in the the interval 
(2^77, z/*) C I~ if jii{Q)At < {2n + l)7r, and with the interval (z/*, {2n + 2)ti) C /+ if 
jii{Q)At > {2n + l)7r, where z/* determines Wmax via p9|) . (This is the vi solution that 
is furthest from {2n + l)vr.) This z^i-value determines a point in the (jir, 7)-plane that 



converges to Pi as Q ^ Qin+i since w(z/i,Q) changes continuously on that subinterval 
as Q ^ Q2n+i- To see that the second of the pair of vi solutions describes a point that 
converges to P2 as Q ^ Q2n+i, we first note that the z/i value must lie between u* and 
(2n + 1)tt and that v* {2n + l)7r as Q ^ Q2n+i- Hence this solution ui (2n + l)7r 
as Q ^ Q2n+i, and the associated point in the (jir-, 7)-plane must approach the line 
jir = —7 as Q ^ Q2n+i since ji,. = 7cos(z/i). Moreover, the value of jir = /(Q) + 7 
converges to f{Q2n+i) +7 as Q ^ Q2n+i; it then follows that this point converges to P2. 
Thus we have shown that this critical curve in the (jir, 7)-plane, parameterized by Q, 
is continuous at any points Q = Q2n+i associated with the parameterization. □ 

The following lemma provides a sufficient condition for a nondegenerate critical 
curve to lie inside of the existence boundary for fixed values of bi, 63, K, and p, i.e. as 
in Fig. infb). Thus, it determines a sufficient condition for stabilization via time-delay 
feedback to fail. 

Lemma 8. Assume that there is a Q-interval Iq = [Qi,Q2] on which ju{Q)^t 7^ 2n7r 
for any n E Z and on which f{Q) > with f{Q) = only at the endpoints Qi and Q2- 
If there exists a value Q E Iq for which D{Q) = 0, where D{Q) is given by (^jj, then 
the traveling wave cannot be (linearly) stabilized for any value of the time delay feedback 
parameter 7. 

Proof. Assume that the conditions of the lemma are met and let Q* be the largest value 
oi Q E Iq for which D{Q) = 0. That the traveling wave cannot be stabilized for 7 > 
follows from part (b) of Lemma El For 7 < we show below that there is a continuous 
critical curve, parameterized by Q G [Q*, Q2], that extends from 7 = jir = to 7 ^ —00 
along the asymptote jir = 7. 

To see that there is a continuous critical curve parameterized by Q € [Q*,Q2], 
we must consider the mapping ()25II2()|1 between ui and Q over this range of Q-values. 
Since D{Q) is continuous, with D{Q) = only at the endpoint Q = Q*, then the sign 
of D{Q) for Q G {Q*,Q2] is the same as the sign of D{Q2). Since WmaxiQ)^ defined 
by (jSni), is positive whenever iii{Q)/S.t 7^ 2n7r, and since /(Q2) = 0, we know that 
D{Q2) < 0. While for Q = Q* there is a unique value ui = u* associated with the 
mapping ()25ll2(jj) . which is given by (jSHI), there are two corresponding values of z/i for 
all other values of Q in the interval since D{Q) < {i.e. for Q G {Q*,Q2))- (If ever 
jii{Q)At = (2n + 1)tt, then one of these z/i values is (2n + l)7r.) All of the values of 
z/i lie in a //-interval = (2n7r,2(n + l)7r) for some n since jii{Q)At 7^ 2m7r for any 
m G Z and any Q G [Qi, ^2]- We may now define two functions i^i{Q) and i^i{Q) using 
the mapping on [Q*, Q2]. We choose I'liQ) to be the solution of ()25ll2(jj) that is furthest 
from the midpoint of Jy {i.e. furthest from Ui = {2n + 1)tt) and i^KQ) to be the solution 
that is closest to ui = {2n + l)7r. We note that the function v^iQ) is discontinuous 
at any Q where jii{Q)At goes through (2n + l)7r due to the jumps in the interval Jj, 
from to 1^7. (In particular, the discontinuity corresponds to being reflected to 
2{2n + 1)71 — u^] see Lemma IHl) In contrast, there are no discontinuities in I'llQ); it 




Figure 8. Schematic of the mappings between vi and Q that parameterize a critical 
curve in the (jir, 7)-parameter plane. The parameterization by Q S [Q*,Q2] is 
determined by the two functions i^iiQ) and Vi{Q), defined in the proof of Lemma|Hl 
where e [2mT,2{n + 1)tt]. The schematic is typical of the degenerate case where 
jii{Q)At — {2n + 1)71 for some Q £ {Q*,Q2), in which case i'^ is discontinuous as 
shown. 

takes on the value (2n + l)7r at any Q value where ju{Q)At = {2n + l)7r. See Fig. |H]for 
an example of functions ViiQ) and v^iQ) in the discontinuous case. Note that the two 
functions approach the unique value i/i = z/* as Q ^ Q* . 

A Q-parameterized critical curve is obtained by inserting the pairs (Q, J-'iiQ)) and 
{Qi ^liQ)) into the parametric equations for ji^ and 7 given by The continuity 

of the critical curve associated with v\{Q) follows directly from the continuity of the 
mapping in this case. The function i'i{Q) may be a discontinuous function at isolated 
points Q where jij(Q)At = (2n + l)7r, as described above. However, by LemmaEl these 
points do not lead to corresponding discontinuities in the critical curve. Thus the portion 
of the critical curve parameterized by (Q, i^l{Q)) is also continuous. The two functions 
together parameterize a continuous critical curve from the origin in the (jir, 7)-plane 
to 7 —00 along 7 = jir. This follows since the two curves meet as Q — > Q*, and 
since {Q2,i^\{Q2)) = (Q2, jij(Q2)At) and {Q2,i^i{Q2)) = (Q2,27rm), where m = n or 
m = (n + 1). Thus the curve parameterized by {Q, i^iiQ)) approaches the origin of the 
(jir, 7)-plane as Q ^ Q2 and the curve parameterized by {Q,i^i{Q)) approaches, as 
Q ^ Q2, the asymptote 7 = jir with 7 ^ — 00. □ 

We now use the results of the previous lemmas to prove our main stability result 
stated in Theorem ^ 

Proof, of Theorem^ The necessary condition for linear stability of the traveling wave 
follows from Lemma |H1 in the case that jii{Q)At 7^ 27m for any n & Z and for any Q 
for which f{Q) > 0. Note, in particular, that the equations ()34IIH5j) in Theorem [T] are 



equivalent to the condition D{Q) = in Lemma |H1 through the definition of WmaxiQ) 
given by (jH^ . 

If jii{Q)At = 27m for some Q = Q2n where f{Q2n) > 0, then Lemma IH states that 
the travehng wave cannot be made hnearly stable (with the prescribed value of p) for 
any value of 7 associated with the time-delay feedback. We now need to show that in 
this case fl34H35|) has a solution. If f{Q2n) = 0, then {Q, ui) = {Q2n, 27m) solves ()34II35|) . 
If f{Q2n) > 0, then we can show there must be a solution to D{Q) = on the interval 
{Q2n,Q2), where Q2 is the zero of f{Q) that is closest to Q2n, i-G. /(Q2) = and 
f{Q) > for all Q G {Q2n,Q2)- (Recall that the equation D{Q) = 0, where D{Q) is 
defined by (jUj), is equivalent to ()34ll35j) .) To see that there is a Q G {Q2n, Q2) for which 
D{Q) = 0, note that D{Q2n) > since Wmax{Q2n) = and D{Q2) < since /(Q2) = 0. 
The result follows from the continuity of D{Q). 

To complete the proof we show that the traveling wave may be stabilized at 7 = — 
if there is no solution of ()34II35|) . This follows since the equations ()34ll35p describe points 
on a critical curve within the accessible region for 7 = — Thus there is no critical 
curve that passes through 7 = — if there is no solution to ()34H35p . In this case, the 
stability boundary must cross the existence line as in Fig. IHl^a), and the traveling wave 
can be stabihzed for some range of 7 values, including J = □ 

The significance of Theorem Q is that it determines an "optimum" value of 7 to use 
in stabilizing the traveling wave; this is the value 7 = ^-^^ To be more precise, any 
set of 7 values that stabilize the traveling wave must include 1 = We note that 

this 7-value coincides with where the stability boundary is furthest from the asymptote 
jir = 7 for each fixed value of juAt 7^ 2mT (see part (1) of Lemma Q). The power of 
Theorem n is that it replaces the "brute force" approach of determining whether there 
are any growing solutions of the linear delay equation (fT^ for each and every value of Q 
for which f{Q) > with the simpler problem of determining whether there is a solution 
to an algebraic equation for Q over the same interval for which f{Q) > 0. 

4. Numerical Linear Stability Results 

Figure ini presents numerical linear stability results in the (|i^|, p)-parameter plane for 
various values of (61,63) in the Benjamin- Feir unstable regime. These diagrams were 
created using the stability criterion of Theorem [T] Specifically, the criterion was apphed 
at points on a grid in the (i^, p)-plane with spacing of 0.0025 in the ii'-direction and 
0.001 in the p-direction. In order to determine whether there is a solution to (jMl - 
l35|) of Theorem 21 we first eliminated the variable u from ()34II35|1 yielding a single 
equation for Q. Care was taken with the trigonometric terms in z/, which needed to 
be inverted so that the u values were in the interval I^, that contained jii{Q)At for 
each value of Q. (See Lemmas EHHl for a description of I,,.) Newton's method was then 
used to determine whether there was a solution for any Q satisfying < f{Q)At < 2. 
Specifically, for each traveling wave tested, the program first checked all values of Q 



between .002 and 20 at an increment of .002, determining the ranges of Q values for 
which f{Q) > 0. (Although we did not prove in all cases that there cannot be an 
instability for Q > 20, this represents an extremely conservative estimate of the largest 
Q value for an instability for the 63, K) values considered; instabilities were generally 
bounded well below Q = 20 for the parameters used in the plots.) Next, on each interval 
for which f{Q)At G [0,2], up to twenty-four Q values were chosen as initial conditions 
for Newton's method, which was used to determine whether there was a solution to the 
stabihty criterion equations ()34IIH5j) or not. 

Note that the black region in Figure Efg) indicates parameter values for which 
f{Q) < for all Q tested (z.e. for Q G (0,20) at intervals of .002). Thus, by LemmaEl 
the traveling wave can be stabilized by spatially-translated feedback alone in this region. 
The lighter shaded regions of Figure 121 indicate parameter values for which the traveling 
wave could be stabilized by a combination of spatial and temporal feedback, at least for 
the "optimal choice" of temporal feedback parameter 7 = The feedback failed to 

linearly stabilize the traveling wave at all other points in the parameter plane. 

Figure ini^h) summarizes our linear stability results in the (61, 63)-parameter plane; 
it indicates both parameters for which we found stable regions in the (i^, p)-plane (c/. 
Figure inta)-(g)), and parameter values, marked by triangles, where either no stable 
traveling waves were found or only very narrow stability regions near the unphysical 
limit in which At 00. From Figure IHl it is evident that the stable regions in the 
(i^', p)-parameter plane are more substantial close to the Benjamin-Feir line. Moreover, 
from a comparison of Figure Efh) with Figure ^ we see that the feedback technique 
seems to be more effective in stabilizing traveling waves in the phase turbulent regime, 
although it also works for some parameters in the amplitude turbulent regime. 

We have performed numerous "brute force" checks of the stability results 
summarized by Figure IHl in order to have confidence in our numerical implementation of 
the stability criterion of Theorem ^ An example of such a check is given in Figure EH 
which shows the results of our direct calculation of the growth rates of solutions of the 
linear delay equation (for j = 1) as a function of Q for two different points in the 
parameter plane associated with Figure IHfg). Figure ITUT a) corresponds to a point just 
inside of the stabihty boundary, while Figure ITOf b) corresponds to a nearby point that 
is outside of the stability boundary. 

We now examine in greater detail some of the features evident in Figure IHl We 
first note that in each of Figures IHfb-g) the leftmost stability tongue extends to p = 
indicating that the associated i^-values can be stabilized with temporal feedback alone. 
We find that the i^-values associated with the left-most stability tongue in FiguresEta)- 
(g) include a value K for which spatially-resonant perturbation wavenumbers Q2 = 2K 
have associated frequencies jii{Q2) that satisfy jii{Q2)At = — vr. The significance of this 
latter relation is that it corresponds to the degenerate situation for which the stability 
boundary in the (ji^, 7)^plane extends all the way to the line 7 = —jir for 7 = — 1/At; 
see Fig. infa) with a = jir, G = 'j. In particular, while perturbations with the resonant 
wavenumber Q = Q2 are not affected by the spatial feedback, their growth rates are, 



= (1.5,1.0) 




Figure 9. (a)-(g) Examples of stability diagrams in the (|i4r|, p)-parameter plane, 
created using the stability criterion of Theorem ^ for the values of (61,63) indicated 
above each plot. The lighter shading indicates regions that can be stabilized by a 
combination of spatial and temporal feedback (with 7 = — 1/ At) . In the darker region 
in (g), the traveling wave can be stabilized by spatial feedback alone, (h) Summary of 
our stability results in the (63, 6i)-parameter plane. The letters refer to the stability 
diagrams shown in plots (a)-(g), which exhibit stable regions in the (A', /9)-parameter 
plane. At points marked by triangles, either no stable traveling waves are found, or 
only very narrow stable regions near the unphysical limit in which At —f 00. 
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Figure 10. Plots of the growth rates associated with solutions of the delay equation 
(PI (for j = 1) as a function of Q for 63) = (2.5, 2.0), p = .007, and (a) K = .2800, 
(b) K = .2875. 



loosely speaking and in the sense just described, maximally reduced by the temporal 
feedback. These observations suggest a possible explanation for why this first stability 
tongue can extend all the way to p = 0, as well as for its positioning along the K-axis. 

Note that the lone stability tongue of Figure Efd) narrows as p increases, eventually 
vanishing (at our resolution in K) once p exceeds 0.04. A similar, though less 
pronounced, narrowing of the left-most stability tongue, with increasing p, is present in 
the other stability diagrams presented in Figure IHl We conjecture that this narrowing is 
related to the observation above that these stability tongues are associated with values 
of K that are unstable in the absence of feedback to spatially resonant perturbation 
wavenumbers Q„ . To understand this claim, expand the real and imaginary parts of the 
eigenvalue ji{Q) of (fTK|l about the resonant Q = Qn = nK (n = 2 above), and focus 
on the specific contribution to ji{Q) that is due to the spatial feedback term. Taylor 
expanding the diagonal factor p(e**^^^ — 1), with Q = Qn + AQ, about AQ = 0, we find 
that the effect of the spatial feedback term is to move the eigenvalue ji in the complex 
plane by an amount p{iAQAx — AQ'^Ax'^/2 + ■•■). For large enough p we thus expect 
the width of the band of unstable wavenumbers about Qn to scale as 1/ ^/p since the 
spatial feedback will decrease the growth rate f{Q) near Qn by pAQ^Ax'^/2; i.e. the 
larger p is, the smaller the unstable band about Qn will be. This in turn leads to an 
0{y/p) contribution to the change in ju over the same band of unstable wavenumbers 
since the spatial feedback shifts ju by pAQAx. From Lemma |3] it follows that if the 
change in juAt over the unstable band exceeds 27i then the traveling wave is necessarily 
unstable. Thus, we expect the stable regions associated with spatially resonant Q must 
vanish for large enough p. 

We now focus on the accumulation of narrow stability tongues that are evident, 
for example, in Figures inte)-(g). We find that the stability tongues become narrower 
and narrower, shifting to larger and larger values of p, as \K\ — > -7===, which is 
where the frequency u ^ and consequently At diverges. Figure fTlT a) shows a blowup 
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Figure 11. (a) A blow-up of a portion of the stability diagram for {bi, 63) — (2.5, 2.0) 
(c/. Fig.EIg)). (b) Corresponding plot cos{jii{Qjnax)^t) vs. K, where Qmax{K) is 
the perturbation wavenumber that gives the largest growth rate fmax for each K. 



of part of Fig. El^g). The rapid alternation between stable and unstable regions with 
small changes in K, occurs as K approaches the singular point where a; — > and 
At — > oo. We trace this rapid variation in stability to the rapid winding of the phase 
jij((5)At, through successive multiples of 27r, when At is large. This claim is explored 
in Figure ITlT b). Specifically, for this figure the value of p is held fixed at p = 0.1, and 
the value of Qmax as a function of K is determined; here Qmax corresponds to the value 
of Q at which the growth rate f{Q) reaches its absolute maximum fmax- Figure ITlT b) 



is a plot of cos{jii{Qmax)'^t) vs. K, which shows that the spacing of the stabihty 
tongues in Figure ITlT a) is the same as the spacing of the maxima of cos{jii{Qmax)^t) 
in Figure ITlT b). On closer examination, we find that the positions of the stability 
tongues in Figure ITlT a) for p = 0.1 are approximately centered around the points for 
which jii{Qmax)^t = (2n + l)iT, for successive integers n. That the locations of the 
stability tongues appear to coincide with the places where ju{Qmax)^t = {2n + l)7r is 
perhaps not surprising. In particular, as already noted, it is at this degenerate value 
that the stability boundary in the (jij, 7) parameter plane extends the furthest from the 
asymptote 7 = ji^, and for this reason it accommodates the largest value of f^ax stably 
(c/. Fig. infa)). (See also Just et al. '^2:^.) Moreover, we know that points for which 
jiiiQmax)^^ = 2n7r cannot be stable by Lemma El so it is expected that the stability 
tongues, if present, would avoid the even multiples of vr. It is interesting to note that the 
spatial length scale for the alternation between stable and unstable regions observed in 
Figure ITlT a) is a direct result of the introduction of time delay into the problem. When 
no time delay is included, the stability of the traveling wave depends only upon /(Q), 
the largest growth rate associated with an eigenvalue of the stability matrix J, but once 
the time delay is added both the real and imaginary parts of ji can influence stability 
in a complicated manner captured by our main stability result. Theorem ^ 

While the presentation of our analysis focused on the case where Ax and At are 
chosen so that Kls.x = |ti;|At = 2tt, it is readily extended to the situation where 
a; At = 2m7r and KAx = 2mT for nonzero integers m and n. FigurefT^presents numerical 
linear stability results obtained for various integers m, n. It is perhaps not surprising 
that the effect of increasing the delay time or the spatial shift is to reduce the extent of 
the stability regions in the (1/^1, p)-plane. Afterall, the larger the value of At, the more 
likely it is that juAt will reach a multiple of 2tt over an unstable band of Q-values. 
Similarly, the larger the value of Ax, the more likely that QAx is a multiple of 2tt over 
an unstable interval of perturbation wavenumbers. However, increasing the spatial shift 
by the factor n does not appear to significantly affect the first stability tongue (compare 
Fig- El (a) and (b) with Fig. IHfe)). This observation is consistent with our association 
of this first stability tongue, which extends all the way to p = with values of Q for 
which QAx is a multiple of 27r; the stabilization achieved over that band of K is due 
to the temporal feedback. Increasing the time delay by the factor m has a much more 
significant impact on this first stability tongue (see Fig. IT^ c) and (d)). It both narrows 
the stability tongues and appears to shift them to smaller \K\. This effect is presumably 
related to the overall increase, with At, in the winding of the phase jii{Qmax)^t over 
an interval of K values, thereby accounting for the shift in the positions of the stability 
tongues (c/. Fig. ITT|). 

5. Conclusions 

In this paper we have examined the effectiveness of the noninvasive feedback control 
scheme, proposed first in the setting of a nonlinear optics problem for stabilizing 
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Figure 12. Examples of stability diagrams for (61,63) = (2.0,1.5) and ojAt = 2rmr, 
KAx = 2n7r, where (m, n) are indicated above each figure. The stability diagram 
associated with the standard case (m, n) = (1, 1) is presented in Fig. |Sfe). 



traveling wave solutions of the one- dimensional complex Ginzburg Landau equation. 
Our approach is based on a linear stability analysis, which involves examining the 
solutions of a linear (complex) delay equation whose coefficients depend on the 
perturbation wavenumber Q. The traveling wave solution is stabilized by the feedback 
if all solutions of the Q-parameterized family of delay equations decay. Our analysis 
leads to a single stability criterion. Theorem Q that allows us to determine stability of 
a traveling wave against perturbations of all wave numbers by determining whether a 
nonlinear algebraic equation for Q possesses a solution or not. Perhaps more significant 
is that our main stability theorem specifies a value of the feedback parameter 7, namely 
7 = — which is guaranteed to work, if stabilization is at all possible. Thus, loosely 
speaking, one of the "control knobs" can be set in advance, thereby minimizing the 
amount of control parameter space that need be scanned in trying to stabilize traveling 
wave solutions. 



In addition to our necessary and sufficient stability criterion of Theorem Q we 
determined simple sufficient conditions under which the control scheme will be ineffective 
and interpreted these in terms of resonance conditions. For instance, we find that if 
there is an unstable wavenumber Q that is an integer multiple of the wavenumber K 
of the targeted solution, then the spatially-translated feedback cannot eliminate this 
instability. Likewise if there is an unstable wavenumber and the frequency associated 
with that perturbation, jii{Q), is an integer multiple of the frequency uj of the underlying 
traveling wave state, then the temporal feedback is ineffective in suppressing the 
instability |22]- A simple consequence of these observations is that it is not possible to 
use the temporal feedback to stabilize the spatially uniform oscillatory pattern at i^' = 0; 
this follows quite generally from the fact that the Benjamin- Feir instability is a longwave 
instability of a translation invariant problem. Specifically, translation symmetry forces 
a neutral mode to exist at Q = and the longwave nature of the instability ensures 
that there are instabilities in a neighborhood of Q = 0. For K = 0, and Q sufficiently 
small, these instabilities are associated with purely real eigenvalues of the linear stability 
matrix J. Such steady modes of instability cannot be suppressed by spatially-translated 
feedback if = 0, and they cannot be eliminated by the temporal feedback either since 
jiiiQ) = if the eigenvalues are purely real. This observation is interesting in light of 
similar results obtained by Harrington and Socolar jTH|; building on results in [23 El] 
they show that the extended time-delay autosynchronization approach to controlling 
traveling waves of the CGLE in one-dimension J7] necessarily fails in two-dimensions 
due to symmetries that force the Floquet multipliers to be purely real. 

For various values of (&i,&3) in the Benjamin-Feir unstable regime, we have used 
our stability criterion to determine numerically where in the (li^l, p)-parameter plane 
the feedback control scheme is effective in stabilizing traveling wave solutions with 
wavenumber K. We find that the stability regions take the form of stability tongues, 
and we offer some insights into the position and spacing with \K\ of these stable regions, 
as well as their extent in the feedback parameter p. While we find some limited regions 
where temporal feedback alone works (p = 0) and another case where spatial feedback 
alone works (7 = 0), most stable regions require a combination of spatial and temporal 
feedback to stabilize the traveling wave. However, a notable feature of the stability 
diagrams in Figure IHl is that the minimal p required to stabilize the waves is often quite 
small. 

We expect Theorem ^ will have applications beyond the example considered here 
of feedback control within the setting of the complex Ginzburg Landau equation. In 
particular, we expect our analysis to apply to other linear stability problems that lead 
to a parameterized family of delay equations of the form ()19p , provided that the growth 
rates jkr{Q), in the absence of temporal feedback, have an absolute maximum as a 
function of Q. Our analysis of the CGLE with feedback was further simplified by the 
fact that we could eliminate one of the delay equations from consideration altogether. 
We did not, however, consider the case where the feedback parameters 7 and p are 
allowed to be complex. 



This paper suggests a number of directions for further investigation. For example, 
it would be interesting to carry out numerical investigations to determine the evolution 
of the system in the parameter regimes where the traveling waves are not stabilized, 
i.e. choose Ax and At appropriate for a traveling wave of wavenumber K that is 
not stabilized by the feedback. Would the system then evolve to a different traveling 
wave state, one for which the feedback does not vanish? More generally, numerical 
investigations, guided by our analysis, could help assess how well the feedback works 
when stabilization is possible, but when it is not possible to prepare the system in a 
neighborhood of the targeted traveling wave pattern. For instance, will the feedback 
work when it is "turned on" after a spatio-temporal chaotic state has developed? 
Another avenue for numerical investigation stems from a potential drawback associated 
with the feedback scheme investigated here - it assumes that the traveling wave 
dispersion relation uj{K) is known. Hence it would be worthwhile to investigate 
numerically the consequences of a "mismatch" between the spatial shift Ax and the 
temporal delay At. This might be particularly interesting in the case of the two- 
dimensional CGLE, in which case the traveling wave pattern can "rotate" in order to 
compensate for any mismatch. 

Another direction for future research is suggested by our choice of the CGLE for 
the detailed analysis presented here. Specifically we refer to the fact that the CGLE 
arises as the universal amplitude equation describing the long spatial and slow temporal 
evolution of homogeneous oscillations in spatially extended isotropic systems. Are these 
the appropriate scales on which to be applying the feedback? If so, how should the 
feedback parameters for the "original problem" be chosen to achieve the optimal value 
of 7 = — 1/At that we find for the "reduced model" given by the CGLE? As a first step 
in addressing these questions it would be worthwhile to carry out an analysis, similar to 
the one we have done here, but on the original equations describing the physical system 
undergoing a Hopf bifurcation, e.g. on a model of chemical reaction-diffusion system in 
the oscillatory regime. 

An attractive feature of the feedback control scheme investigated here is its 
noninvasive nature, i.e. the feedback vanishes when the targeted state is reached. This 
design feature exploits the underlying spatial and spatio-temporal symmetry properties 
of the desired pattern in a very natural way, suggesting that this feedback approach 
may be readily generalized to more complicated patterns. For instance, the original 
problem investigated by Lu, Yu, and Harrison |2] involved oscillatory stripe patterns 
in a two-dimensional isotropic system. They found that they could target a particular 
orientation of the stabilized stripe pattern by an appropriate rotation of the spatial 
translation applied in the feedback control. Thus it would be of interest both to extend 
our linear stability analysis to the control of traveling plane wave solutions of the two- 
dimensional CGLE, and to consider the feedback stabilization of more complicated 
traveling wave patterns. 
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